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Abstract 

We demonstrate that the efectric field enhancement that occurs in a gap between two closely 
spaced nanostructures, such as metallic nanoparticles, is the result of a transverse electromagnetic 
waveguide mode. We derive an explicit semianalytic equation for the enhancement as a function 
of gap size, which we show has a universal qualitative behavior in that it applies irrespective of 
the material or geometry of the nanostructures and even in the presence of surface plasmons. 
Examples of perfect electrically conducting and Ag thin-wire antennas and a dimer of Ag spheres 
are presented and discussed. 
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FIG. 1: Schematic diagram of a thin-wire antenna. The parameters shown are discussed in the 
text. 

Structures that generate large electric field enhancements relative to the incident field, 
hereon referred to as |E| 2 enhancements, have recently received a great deal of attention 
[THI0]. This is because such enhancements are central to a number of physical processes, 
including such surface-enhanced spectroscopy techniques as surface-enhanced Raman scat- 
tering (SERS), second harmonic generation, and enhanced absorption and fluorescence pp. 
Often the enhanced field is generated in the small crevices of a roughened metal surface 
[2] or at the junctions of closely spaced nanoparticles [3TTTU]. Herein we focus on the latter 
structures, and while much is known about the |E| 2 enhancements in them, there is still 
confusion over fundamental principles. In particular, the functional dependence on gap size 
[3 El El [9] , arguably the most basic and important aspect, has not been quantitatively de- 
termined and the underlying physical principles which determine it are not entirely known. 
It is the purpose of this Letter to resolve this issue through finite element method (FEM) 
calculations [TT] and an analytical theory developed for the transmission of light through an 
isolated slit in a metal film [T2] . 

For two closely spaced nanostructures, the |E| 2 enhancements in the resulting gap can in 
principle be explained using antenna theory [131 EH] > where the open-circuit voltage across 
the gap is responsible, and thus the systems are often classified as such [21 H]. We therefore 
begin by considering a two-dimensional (2D) antenna as shown schematically in Fig. [T] [The 
extension to three dimensions (3D) will be discussed below.] Two metal wires M with widths 
w are separated by a distance a and the entire structure spans a length of h. The structure 
is illuminated from below at normal incidence by a plane wave with wavelength A, and we 
wish to determine how the |E| 2 enhancement at the center of the gap depends on a. It 
is important to realize that for a real metal and distances less than approximately 1 nm, 
nonlocal dielectric effects will become important [15J. Our quantitative analysis herein will 
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FIG. 2: |E| 2 enhancements as a function of gap size for PEC thin- wire antennas with parameters 
given in the text. Solid lines are used to connect the actual data points (symbols). 

therefore be for a > 1 nm, but in most cases we will include smaller distances to highlight 
qualitative features. Antenna theory for a perfect electrically conducting (PEC) thin-wire 
antenna (w <C h) assumes that the incident electric field E generates an alternating current 
along x, which results in an induced voltage V across h. If h ~ n\/2, where n is an integer, 
the antenna resonates and V ~ —Eoh, where Eq is the amplitude of the incident field [IB]. 
As a result, the open-circuit voltage in the gap should produce a uniform |E| 2 enhancement 
of |E|7|E | 2 = \V/a\ 2 /\E \ 2 « h 2 /a 2 . 

In order to test the above analysis, we rigorously determined the |E| 2 enhancements via 
FEM calculations [TT] for h = 250 and 500 nm thin-wire antennas (n = 1 and 2, respectively) 
with w = 5 nm at A = 500 nm for gap sizes of a = 0.125 to 10 nm. To characterize the a 
dependence, we can assume that |E| 2 /|Eo| 2 is proportional 1/oP and plot the results on a 
log-log scale to determine p; Fig. [2j For both antennas it is found that p « 1.2 for a > 1 
nm, and it is even less for smaller a, which is much lower than the above antenna theory 
prediction of p = 2. An alternative way to describe these systems and the |E| 2 enhancements 
that they exhibit is therefore needed. 

The system in Fig. [l] can be greatly simplified by taking h — > oo, which below we will 
show does not significantly affect the behavior of the electric field E in the gap. Maxwell's 
equations can be solved analytically for such a system (if the metal is a PEC), which is an 
isolated slit in a metal film, by appropriately expanding the transverse component of the 
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field (the y component of the magnetic field, in this case) above and below the film and inside 
of the gap in terms of known functions and applying boundary conditions at the interfaces. 
While the full solution for this problem has been implicitly worked out in terms of a system 
of linear equations [12J, we demonstrate that under a few reasonable approximations it is 
possible to obtain a tractable semianalytical form for |E| 2 /|Eo| 2 . 

Inside the gap E can be defined entirely in terms of its x component E x (below we will 
show that the z component E z is zero), which can be expanded as a superposition of forward 
and backwards propagating (and evanescent) waveguide modes m, 

oo „ 

E x (x, z) = T^ ( A me^ mZ - B m e'^ z ) <j> m (x) , (1) 

m=0 u 

1 /2 

where A m and B m are the respective modal amplitudes, ko = 2n/\, (3 m = [k$ — (imc/a) ] , 
and (j) m (x) = (2/a 1 / 2 ) cos [inn (x + a/2) /a] is the solution to the Helmholtz equation subject 
to PEC boundary conditions on the gap sides at x = ±a/2. A m and B m can be found by 
ensuring the continuity of Eq. ([!]) at the input (I) and output (O) surfaces of the gap, 
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where E^ and E^ are intensities of scattering events that take place at each of the surfaces, 
which can be determined by solving a set of linear equations (see below) [12] . 

Inside a gap that is small relative to the incident wavelength (/c a <C 1), the only waveg- 
uide mode that can exist is the m = transverse electromagnetic (TEM) one [T^. In this 
case, there is no x dependence in E x , E z is zero, (3 = k , and the linear equations defining 
Eq and Eq simplify considerably, 

£j = /o 5 (4) 

(/o + ^oo) 2 - {gtf 

K = -j^rEi (5) 

Jo + 9oo 

where Jo = AEqcl 1 ! 2 (for normal incident light) is the overlap amplitude of the incident field 
with the TEM mode, f = % cot (k w) is the admittance amplitude, Qq = i csc (k w) is the 
coupling amplitude between I and O, and 

4 W 2 ra/2 

9oo = - / dx dx' G (x, x') (6) 

a J -a/2 J -a/2 



is the amplitude of the TEM mode's self-interaction, where G (x, x') = 
(ko/2)H^ {ko\x — x'\) is the 2D vacuum Green's function with Hq 1 ^ being a Hankel 
function of the first kind. The maximum value that k \x — x'\ can take is k a. Since 
koa < 1, we can make the small-argument approximation in Hq and perform the integral 
in Eq. ^ analytically, 

goo = 2k a(l + il) (7) 

where / = (2/tt) [In (koa/2) + 7 — 3/2] with 7 being Euler's constant. 

For a gap with a small width relative to the incident wavelength (kow <C 1) we can use 
the small-angle approximation in Eqs. (|2j) and ^ to greatly simplify Eq. ([I]), 
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Note that the approximation g 00 <C 2f was also used to get Eq. ([8]), which is valid consid 
ering that the leading terms are koa and 1/koW, respectively. 

At the center of the gap (z — w/2) Eq. (J8l) simplifies even further, 



E x (w/2) = 8E -^— (9) 
25-00 + u 

where u = fo — (go) 2 / fo — — ^tan (kow) ~ —ikow. Equation ^ shows that at the center of 
the gap E is inversely proportional to the interference between two terms, the self-interaction 
term goo , which depends only on koa, and a term u representing the interference between 
surfaces I and O, which depends only on k w. 

Using the explicit expressions for g Q and u in Eq. ^ and calculating |E| 2 /|E | 2 gives 

64 

|E|7|E | 2 = \E x (w/2)\ 2 /\E \ 2 = 2 . . (10) 

11/11 ' V 1 n/l ' 16 (1 + 1 2 ) (k a) + 81 (k a) (k w) + (k w) V ' 

For gaps with a not very small relative to w (greater than approximately 1 nm for 
the thin-wire antennas discussed above), E depends primarily on g o and |E| 2 /|E | 2 ~ 
64/ [16 (1 + I 2 ) (koa) 2 ^ . Because the enhancement is of the form l// 2 a 2 , we would expect a 
l/a p fit to not work at all. However, over a couple order of magnitude range of a values 
I ~ In (k Q a/2) is well-approximated by A (k a/2) b , where A and b are constants that can be 
determined by demanding that this equality and a corresponding one involving its derivative 
be satisfied at some value of koa/2. With these constraints one finds that b = 1/ln (koa/2), 
which varies slowly with koa/2. For A = 500 nm and 1 < a < 10 nm, b ~ —0.25 and there- 
fore l/l 2 a 2 ~ 1/a 1,5 , which is consistent with p w 1.2 found for the thin-wire antennas in Fig. 



[2| For small gaps both a 2 and l 2 a 2 go to zero, which means that E depends primarily on u 
and |E| 2 /|E | 2 ~ 64/ (fc w) 2 . Therefore, for a « w we expect a turnover to a weaker l/a p 
dependence, which is also consistent with the results in Fig. \2\ Actual |E| 2 enhancements 



calculated using Eq. (10) for a = 0.125 to 10 nm, w = 5 nm, and A = 500 nm are shown in 
Fig. [2] as well, and are in agreement with these remarks. 

The strong agreement between the modal and thin- wire antenna results in Fig. [2] suggests 
that the |E| 2 enhancements in both cases arise from the same effect, a TEM waveguide mode. 
The former does show a slightly stronger l/a p dependence, but this can be understood as 
follows. Recall that the amplitude for coupling incident light into this mode is proportional 
to their overlap; see Eq. Q. In a finite structure, such as a thin- wire antenna, some of 
the impinging incident light can be effectively lost via scattering, leading to a less efficient 
coupling into the TEM mode and a weaker l/a p dependence. Quantitatively, the amount of 
scattering is given by the scattering efficiency Q sc (the ratio of the scattering cross section 
to the geometric one). In the modal results Q sc is naturally 0, since the geometric cross 
section is infinite. For a finite structure, however, Q sc > 0. For example, the resonant 
antennas h = 500 and 250 nm have similar Q sc values of 1.953 and 1.915, respectively, at 
a = 2 nm. For an off-resonance condition we expect less scattering, and in fact this is what 
is numerically found for h = 175 nm. In this case, Q sc = 0.931 nm (at a = 2 nm) and the 
\ja v dependence is indeed stronger than for the two resonant antennas. 

Thusfar we have considered 2D systems. While such structures are experimentally realiz- 
able, in most cases (e.g., typical SERS substrates) structures are 3D in character. Nonethe- 



less, our analysis remains valid and Eq. (10) should still apply (qualitatively, at least). This 
is because the TEM waveguide mode suggested as responsible for the |E| 2 enhancements 
also exists in 3D. The only requirement to support a propagating electromagnetic wave is 
that an oscillating potential difference be established between the walls supporting the wave 
(e.g., the sides of a gap) [T7j. In addition to 3D, actual structures are comprised of a real 
metal, typically Ag or Au at optical frequencies due to possible increases in |E| 2 enhance- 
ments via surface plasmon (SP) excitations [7J [TO]. Based on the discussion above about 
Qsc, it is reasonable to suspect that the dominant effect of this is that there may be wave- 
length dependent modulations to the 1/aP dependence, due to corresponding dependencies 
in the absorption efficiency (Q a b s ) and Q sc [IS]. However, the overall trends and underlying 
physical principles should remain the same. 
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FIG. 3: |E| 2 enhancements as a function of gap size for a 3D Ag thin-wire antenna and a dimer 
of Ag spheres with parameters given in the text. Solid lines are used to connect the actual data 
points (symbols). 

As a first example of the applicability of our analysis to real 3D structures, FEM was used 
to calculate the l/a p dependence at the center of a 3D Ag thin-wire antenna with h = 250 
nm and w = 5 nm (both in and out of the plane of Fig. [I]) at A = 500 nm; Fig. [JJ It is found 
that p is again approximately 1.2, which is nearly equal to the analogous 2D PEC thin-wire 
antenna results in Fig. [2j Further similarity comes from the behavior as a decreases, where 
the l/a p dependence again becomes weaker, as can be inferred from the curvature of the 
actual data relative to the linear fit. 

It is possible to verify the existence of a TEM waveguide mode by looking at profiles 
of E inside the gap, as for normal incident linearly-polarized light this mode has the same 
polarization. Figure [4] shows the fields inside the gap of the antenna discussed above for 
a = 2 nm, and it can be seen that this is indeed the case, as |E| 2 \E x \ 2 . (Note that 
\E y \ 2 and \E Z \ 2 are both less than 1 on the scale in Fig. |4j) In fact, it has recently been 
demonstrated experimentally that electromagnetic fields inside the gaps of nanostructures 
are linearly polarized, even in more complex ones than discussed here |19j . 

As a further and final example of the applicability of our analysis to real 3D structures, 
|E| 2 enhancements as a function of gap size were calculated for a dimer of 250 nm diameter 
Ag spheres at A = 633 nm (a popular type of experimental system and common laser 
wavelength [6]); Fig. [3] The l/a p dependence in this case is found to be characterized by 
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FIG. 4: Intensities of |E| 2 and the incident component of E, l-E^I 2 , in the gaps of a 3D Ag thin-wire 
antenna (top) and Ag sphere dimer (bottom) with parameters discussed in the text. The inset 
white scale bars correspond to 1 nm. Note that the intensity values have been rescaled relative to 
Fig. [3] to fit clearly on the same scale, and the fields inside of the structures have been set to 0. 



p w 1.3, which is nearly equal to the modal result of p ~ 1.4, and slightly greater than the 
somewhat analogous thin-wire antenna result in Fig. [3j The stronger l/a p dependence is 
related to the fact that this structure is more efficient for capturing light [H], as indicated 
by Qabs — 0.200 as opposed to 0.087 for the thin-wire antenna, for example. Field profiles 
inside the gap again indicate the presence of a TEM waveguide mode; Fig. |4| Looking closely 
at Fig. [3] reveals that there is a much less strong turnover to a weaker \jdP dependence for 
smaller a than was seen for any of the other structures. Such behavior is understandable 
considering that w is effectively zero in this case (there is only a single point of minimum 
approach), which can lead to |E| 2 being unbounded as a — > [lim^o (1 + I 2 ) (k a) 2 = 0; see 
Eq. (10)]. It is quite remarkable that the simple analysis derived for a 2D PEC film with an 
isolated slit is so accurate when applied to full 3D structures of other geometries, and even 
in the presence of SPs. 
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